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Abstract 

We study the conformal field theory of a free massless scalar field living on the half line 
with interactions introduced via a periodic potential at the boundary. An SU{2) current 
algebra underlies this system and the interacting boundary state is given by a global SU(2) 
rotation of the left-moving fields in the zero-potential (Neumann) boundary state. As the 
potential strength varies from zero to infinity, the boundary state interpolates between the 
Neumann and the Dirichlet values. The full S-matrix for scattering from the boundary, 
with arbitrary particle production, is explicitly computed. To maintain unitarity, it is 
necessary to attribute a hidden discrete "soliton" degree of freedom to the boundary. The 
same unitarity puzzle occurs in the Kondo problem, and we anticipate a similar solution. 
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1. Introduction 

Conformal field theory can be denned on manifolds with boundaries, provided that 
appropriate boundary conditions are imposed |IJ. The Dirichlet and Neumann boundary 
conditions on scalar worldsheet fields are familiar, if trivial, examples. Non-trivial examples 
can arise from the introduction of interactions localized at the boundary. In fact, many 
important problems in condensed matter and particle physics, including open string theory 
P,[M, monopole catalysis ||, the Kondo model 0, dissipative quantum mechanics 
and junctions in quantum wires can be described this way. 

In this paper we will present an exact solution of the simplest of all such theories: 
a single massless scalar field interacting through a sinusoidal potential localized on the 
boundary. The period of the potential is taken to be such that perturbatively it is a 
marginal operator. If the potential strength is zero, a free (Neumann) boundary condition 
is imposed at the origin, while infinite potential strength leads to a fixed (Dirichlet) condi- 
tion. Both conditions are trivially consistent with conformal invariance and the question 
being asked here is whether there are conformal boundary conditions which interpolate 
between them, i. e. whether our theory is conformally invariant for arbitrary potential 
strength. Some authors have given partial evidence for the existence of this interpolating 



conformal field theory [JTy, JTJ , but a complete construction of the corresponding boundary 
dynamics has been lacking. Significant steps toward this goal were taken in a previous 
paper by two of us where the partition function and some S-matrix elements in the 
interpolating theory were calculated and found to be consistent with conformal invariance. 
In this paper we present a new method that allows us to obtain completely explicit (and 
very simple) expressions for all the dynamical quantities in the theory: partition functions, 
boundary states, S-matrix elements and so on. 

The scalar field may be compactified at any radius consistent with the period of the 
potential, and one finds that all such radii are integer multiples of the self-dual radius, 
R 8 d — \/2. The key to our calculation is the full exploitation of the SU(2) symmetry 
associated with the boundary interaction (at infinite radius this symmetry is somewhat 
hidden, while at the self-dual radius it is fully explicit). One very interesting feature of this 
SU{2) is that it leads to the existence of "soliton" sectors (basically the SU{2) partners of 
the usual massless boson Fock space states) in the scattering problem. These sectors are 
not visible in naive perturbation theory and, unless they are accounted for, the S-matrix 
is non-unitary! A similar problem afflicts the conformal field theory describing the Kondo 
model, and we expect a similar resolution. 
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2. Exact Boundary State by SU(2) Methods 

2.1. Setting Up The Perturbation Expansion 

Our starting point is the Lagrangian for a free field on the segment < a < I, with a 
boundary interaction at a = 0, 

L=^f da{d,Xf - \{ge iX ^'^ + ge~ iX ^) . (2.1) 

In order to avoid infrared complications we impose a Dirichlet boundary condition at a = I 
( X\ a= i = 0), while the boundary condition at a — is dynamical: 

--L^ + ige iX /^ - ige- iX '^ = . (2.2) 
V2tt da 

The potential strength g is taken complex so that we can vary the location of the poten- 
tial minimum. By varying \g\ between zero and infinity, we effectively interpolate between 
Neumann and Dirichlet boundary conditions at the origin. With the particular period cho- 
sen above, the boundary potential has perturbative scaling dimension one (as a boundary 
operator) and the theory is scale and conformal invariant to first order in an expansion in 
the potential strength. We will show that it is exactly conformal invariant for any value of 
the potential strength and will construct conformal boundary condition that interpolates 
between the limiting Neumann and Dirichlet cases. 

Our first task is to compute the functional integral, which, as is well known |13| , can be 
regarded either as an open string partition function Z BD = tr(e~ TH ) or as the amplitude 
for free closed string propagation between two boundary states Z BD = (B\e~ l ( Lo+LoS> \D) 
where \D) is the Dirichlet boundary state and \B) is the boundary state induced by the 
interaction. Expressing \B) as the Neumann boundary state, \N), acted on by the potential 
term in the path integral gives the more explicit expression 

Z BD — (N\e~ J 9 'e~ ( 0+ °'\D) (2.3) 

In what follows we study the closed string representation, returning to a discussion of open 



string physics in section 4. The expansion of (|2.3|) in powers of the potential gives 
z bd = Z BD J2^(f[f dU (-£ e «X*<>°)/V5 _ l e -iX( ti ,o)/V2^ 

o n ' \i=l^° 



(2.4) 

/ 2\ — 1/24 00 
y BD _ Ig ) \^ i i\n2n 2 



n=— oo 
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where q = e~ 27rl / T and f(q) = n^Li(l — ? n )- The brackets are to be evaluated using the 
free-field propagator appropriate to the boundary conditions of the free problem (periodic 
in time, Neumann and Dirichlet at the two spatial boundaries): 

(*(*!, 0)X{t 2 , 0)) a=0 = -2 log e l^-J ]2i ^ (2-5) 

where r = l/T. 

In it was argued that Z BD must have the general form 



:BD{3 ' s) = vm^W) \ 1 + W 2 cos + " A(9 - »1 ) (2 - 6) 



n=l 



and explicit computation up to third order in perturbation theory gave A(g,g) = 7r(g + 

3 .. 

g)/2 + fg((7 3 + g 3 — 3g 2 g — 3gg 2 ) + . . .. In this paper we rederive eq. (|2.6|) and show how 
to do the perturbation calculation to all orders, obtaining simple analytic results in the 
process. It is convenient to rewrite (|2.6|) in the identically equivalent form 



tBD i 



i=o,i,i, 



sin(2j + l)6>/2 ((g 2 ) 3 -(g 2 ) (3 " +1) ) 

sin^/2 v / 2(g 2 ) 1 / 24 /(g 2 ) 

1 

V2 ; 



j=0,i,l,--- 



(2.7) 



where 6 = 2A(g,g) + tt, Xj U ^ 2 \@) is the spin-j SU(2) character, and xj^iq) = 
f(q)~ 1 q~ 1 ^ 24 (q : ' — q^ +1 ^ ) is a "discrete state" Virasoro character (more about this 
shortly). What is the origin of the SU{2) symmetry manifested in ( |2.7| )? In the con- 
formal field theory of a single free boson on the open line (with no boundary condition), 
the primary fields are: e ipX for any value of p (weight p 2 /2) plus a discrete set of extra 
primary fields that appear at momenta p = nj\f2 [|14| , [T5[| . These discrete states are orga- 



nized in SU(2) multiplets and can be denoted *&j j z where j z is related to momentum by 
j z = p/V2. They have conformal weight h = j 2 , and their Virasoro characters are the 
Xj r {q) defined above. It is clear from (|2.7|) , and from the calculations in [12], that the 
discrete states saturate the dynamics of our boundary problem. The reason is that the 
Neumann boundary state has momentum zero, while the period of the boundary potential 
is such that it injects momenta which are integral multiples of 1/ \[2\ precisely the values 
carried by the discrete states. Other momenta are possible, but they are not excited by 
the boundary interaction. 
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The chiral SU(2) generators which transform the holomorphic discrete states of a given 
j among themselves are 

j± = I ^ e ±iV2X(z) j3 = I ^J_ ldX ( z ) (2.8) 



2ni J 2m 

(there is a corresponding set acting on the anti- holomorphic states). They commute with 
all Virasoro generators because they are contour integrals of weight-one fields. Since X(z) 
is the holomorphic part of the field X, are well-defined only when acting on states of 
momentum p = n/ y/2, which are in fact the only momenta that appear in (2~77) . That is 
why even the R = oo theory has SU(2) symmetry. 

2.2. Exact Calculations at the Self-dual Radius 

The above remarks show that, as far as calculations of the partition function Z BD 
are concerned, no information is lost if we compactify the boson at the self-dual (SU(2)) 
radius, setting X ~ X + 2n\[2. The SU(2) momenta n/V2 are then the only ones allowed 
(but we have to treat left and right momenta as independent and we have winding states 
with pl 7^ Pr). We will now show that in this compactified theory, the integrals needed to 
evaluate the perturbation expansion of the partition function can be converted into global 
SU{2) raising and lowering operators so that the partition function can be calculated by 
algebraic methods! Later on we will show how to extend our simple SU(2) algebraic 
manipulations to the infinite radius problem. 

Our goal is to determine the exact closed string state \B)su(2) corresponding to the 
insertion of a dynamical boundary. As a first step, we recall that the state \N) SU , 2 \ which 
implements the Neumann boundary condition is 

oo 

i«W = E i^Jvf^ (2 ' 9) 

n= — oo v n * 

The normalization factor 2~^ can be extracted by viewing the closed string partition 
function as the modular transform of the open one. It is a measure of the zero temperature 



entropy of the Neumann boundary condition [16]. As shown in ref. [12|, it is helpful to 
represent ( |2.9| )as a sum over Virasoro modules: For each primary field (p we can build a 
reparametrization invariant boundary state \<p)) = Yl n \4>-> n ) 10; ^) where the sum is over 
all the descendants of (p in the module |PT |. At the self-dual radius, the only primary fields 
are the "discrete states" $ Jim $j jm /, and we denote the corresponding reparametrization 



invariant boundary states by \j,m,m')). By examining the left and right momenta of 
states contributing to ( |2.9|) , we see that only the modules with ml = —m should enter. We 
also note that \N) su / 2 -\ is a sum over all possible normalized (to 2 - ^) oscillator states and 
after we change our basis to Virasoro primaries and their descendants, every contributing 
state should again be normalized to 2 - ^. Therefore, we may choose a phase convention 
where every contributing module enters with coefficient and rewrite the Neumann 
state as 

Now consider the oscillator basis version of the Dirichlet boundary state 0: 

oo 

71=— OO v v ^ 

By the same argument as before, this can be written as a sum of \j, m, m)) modules (which 
don't contribute to the Neumann state except for m = 0). Because of the minus signs in 
( |2.11| ), there will be relative phases but, with a little thought, one can see that phase 
conventions for the \j, m, m)) can be chosen such that 

\ D )su(2) =2-^e- iwj \j,m,m}) = ^ Jl Wsup) ■ ( 2 - 12 ) 

The convention-independent fact that makes this possible is that the \ j, 0, 0)) modules enter 
\D)su(2) an< ^ \N)su(2) w ith relative phase (— l)- 7 , a fact already noted in [T^]. Remarkably, 



the Dirichlet state is just an SU{2) l rotation of the Neumann boundary state! 

At the self-dual radius, we can divide the field X into left- and right-moving parts 

X(t,a) = x° + 2np(t + ia) + X osc (2.13) 

X(t,a) = x° + 2np(t - ia) + X osc (2.14) 

where X osc is the part that contains the oscillators, and the left and right momenta and 
zero modes are independent. The Neumann boundary state (|2.9| ) reflects right-moving into 
left- moving fields as follows 0: 

x°\N) =x°\N) ; p\N) = -p\N) ; a n \N) = a- n \N) ; X(t,a)\N) = X(t, -a)\N) . 

(2.15) 



5 



These properties will enable us to rewrite the action of the boundary perturbation in terms 
of SU(2) generators. 

We demonstrate this by using commutativity of right- and left-movers to rearrange 
the typical product of normal ordered exponentials appearing in an O(g n g ) term in the 
expansion of (|2.4j) as follows: 



n n ^ n 

(N\ Y[ ■ e i{xiti ' 0)+x{u ' 0))/ ^ : = (N\ JJ : e iX ^ )/^ : JJ : e «"(*i,o)/v^ . ( 2 .16) 

i=l i=l i=l 



We then use ([2.15 ) to convert the leftmost right-moving exponential into a left-moving one 



and then commute it to the right of the remaining right-moving exponentials. When all 
the right-movers have been so eliminated, we have 

(N\ : e iX (*»>°)/v^ • • e iX(t!,0)/V2 .. e iX(t u 0)/V2 . . e iX(t n ,0)/V2 . 

The central pair of terms can be combined into the dimension-one current C : e l ^ x ( t ^^) ■ 
where C is a normal-ordering constant that will be absorbed in g. Remember also that 
the the symbols X and X are now being used to denote the two chiral components of the 
field. The integral over the time argument converts the current to the global SU(2) raising 
operator: 

J dte^ x ^ = j dze^W = 2iriJ+, z = e^ u / T (2.17) 

(no Jacobian is needed when changing variables from t to z because we are integrating 
a weight one field). We can freely move J + to the left because the intervening q %x I^ 
operators correspond to (j,m) = (i, |) states and commute with the J + raising operator. 
This maneuver can be repeated until all the potential insertions have been converted to 
577(2) raising operators. The result for the O(g n g ) term is 



V2 n\ *— ' 7 *— ' 

j',m' j,m 



The right-moving states are not affected by the J + operators and their inner product sets 
j = j' and m = m' . The sum over descendants in the Ishibashi states produces a factor of 
xj ir (q 2 ) and all that is left is the SU(2) matrix element (j, | \(J+) n \j, -f ). In the end, 
we will not even need to calculate this matrix element, because the sum over n generates 
something even simpler. 
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But first we have to analyze what happens to the general g n g m term in the expansion. 
In a general term, the operators will encounter e ±lX /^ operators with which they do 
not commute and we will have to keep track of SU (2) commutator terms as we reduce the 
amplitude. Consider the example of an O(gg) term, starting with the stage where the first 
pair of j = \ fields have been combined into a raising operator density: 

dtidt 2 . e - iX (t 2 ,o)/V2 .. e iV2X(t u o) .. e -iX(t 2 ,o)/V2 . = 

(2 " )2 rat (2 - 19) 

(N\J+J- + / g(iV| : e «(*».0)/^ :: e -«(t 2) o)/V2 . 

The last term comes from the commutator associated with moving J + past e ~ lX /^, The 
coincident-point product of e lX l^ with e - tX /^ evaluates to a divergent constant which 
can be absorbed as a constant shift in the interaction potential. This has no effect on the 
physics and can be dropped. When we try to extend this sort of argument to higher order, 
we find another type of divergent commutator which can this time be absorbed as a finite 
renormalization of the coupling strength"!". This procedure can be generalized to all orders 
but we will defer the details of the argument to Appendix A. There we will show explicitly 
how to regularize the theta function integrals so that all commutator terms generated by 
moving J ± operators past j = \ fields can be discarded. 

Given such a regulator scheme, any term of order g n g m reduces to the matrix element 
of a product of n raising and m lowering operators in some particular ordering. The sum 
over all terms with n + m = N is easily seen to give 

Z n = (mr^(N\ sui2) (-gJ+-gJ-rq L o+^\ D)su{2) = (2 .20) 
The final sum over n can be done explicitly, yielding 



? 1 


E 


o,|,i,... 


\m=-j 


? 1 

o,|,i,... 


, E 0', 

\m=-j 



v / 2 / 



(2.21) 



^ To be a bit more specific: With a standard short distance cutoff e, all divergences can be 
absorbed by choosing the "bare" coupling constant to be g/e. The procedure we are outlining 
amounts to choosing the "bare" coupling to be a power series ^(1 + ci|g| 2 + • • -)/e. This doesn't 
change the physics of the theory, but does change the precise meaning of the parameter g. 
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where the second line is obtained using the SU(2) relation e l7rJl \j, m) = e i7rJ |j, —m). 
Since the sum over m gives an SU(2) character, it is easy to see that (|2.21|) reduces 



to fl2/Fp with A defined by 



/ « / n n (q + q) sin(7r|g|) 
sm(A(g,g)) = Ky y) yi> . (2.22) 
^ \9\ 

This expression has a curious implication about the limit of infinite potential strength. We 
expect that when \g\ — > oo the boundary state \B) turns into a sum over Dirichlet states 
with the field sitting at the minima of the potential (whose locations are in turn set by the 
phase of g). This means that, in the limit of infinite potential strength, A should approach 
the phase of the complex coupling g. According to ( |2.22p , this happens at \g\ = 1/2 rather 
than at \g\ — > oo. This is just a finite renormalization effect: As we remarked earlier, 
our coupling is related to the usual one by a coupling constant redefinition of the form 
9' = 9 f {99 s , I • This can have the effect of mapping infinite coupling strength to \g\ = 1/2 
if / has a singularity there. 

The most important fact about ( |2.21| ) is that, comparing with ( |2.9|) and (|2.12| ), we 
can read off the exact interacting boundary state: 

\B)sum = e^ J++ v J ^\N) su{2) = e ieaja \N) su(2) , (2.23) 

— * 

where 6 = 27r[Re(g), —Im(g), 0). The net effect of the interactions is to carry out a global 
SU{2) rotation, through angle 0, of the left-movers with respect to the right-movers in 
the original Neumann state! Since SU(2) has three generators, it may be surprising that 
we don't find a three-parameter family of boundary states. Actually, there is a somewhat 
more general theory with essentially the same boundary dynamics: We may add to ( |2 . 1| ) 
a term of the form aX(a = 0) without destroying the essential features of our analysis 
and the resulting SU{2) rotation depends on three parameters, Re(g), Im(g) and a. For 
simplicity we set a = in the rest of the paper. 

\B) imposes reflection boundary conditions on a) which are global rotations of the 
conditions ( |2.15| ) imposed by \N). We will discuss this in more detail in section 3. The new 
state obeys, as it should, the reparametrization invariance conditions (L n — L— n )\B) = 
because \N) obeys them and the SU(2) generators commute with the Virasoro generators. 
This is a useful consistency check of our regularization procedure. 



' In fact, the regularization used in [|12| is the standard one and the coupling constant g' used 
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there is related to ours by g' = g(l + j^gg + 
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2. 3. Extension to Other Radii and Boundary Conditions 

Now let us consider other compactification radii for the field X. The only radii 
consistent with the marginal potential in (|2.1|) are integer multiples of the self-dual radius. 
We can extend our techniques to solve the theory at any of these radii but the results are 
particularly simple in the R — oo case. Let us consider calculating the partition function 
Z BD at infinite radius. At this radius the Neumann and Dirichlet boundary states have 
the form 

1^=00 = 2"^ a " a "|0> = 2-* Yl b'»M», ( 2 - 24 ) 

j',m=0 

\D) R=oc = 2-Ve- ^ a " a " J dp\p) R \p} L . (2.25) 

The Neumann state carries momenta Pl = PR = 0. Since the interaction term changes 
left and right momentum equally and in multiples of l/y/2, the only states from \D) R=oo 
that will contribute are the ones that have the momenta Pl = PR = n/y/2. But these 
are precisely the momenta that appear in the Dirichlet state at the 577(2) radius, so 
we can, without any error, make the replacement \D) R=OQ — ► \D) su , 2 y We can also 
replace \N) R=QO by \^)su(2) m ^ ne ma trix element since the difference between the two 
is a collection of terms with momenta pl = ~Pr = n/y/2 with n ^ 0, none of which 
have any overlap with the Dirichlet state (which has Pl = Pr)- The net result is that 
Zr=oo = ^su(2)- -ft ^ S eas y to extend this argument to any allowed radius and find that 

yBD _ yBD 
n R ~ Zj SU(2)' 

The boundary state itself is a bit more interesting. The boundary state at the 577(2) 
radius ( |2.23| ) has contributions from all 577(2) values of pl and pr independently. At 
multiples of the SU (2) radius only a sublattice of the SU (2) pl and pn values are allowed 
and one gets the boundary state by projecting ( |2.23|) onto the allowed momenta. At 
R = oo, the condition is pl — Pr and we have 

\ B )r=°o =2-^J2J2\j,n,n))((j : n : n\e ieaja \j,m,-m)) 

j m,n 

(2.26) 

=2_i E E V L-m\j^M) 

j m=-j 

where T> 3 m _ m is the rotation matrix element 

< DL,-m = {jM^ {9J++ " J - ) \j, ~m) 
9 



This expression is identical to the R = oo boundary state found in [ff2[. We do not 



have contributions from the continuum states because \B) is generated from a Neumann 
boundary state with an operator that changes the momentum in discrete steps and thus 
only the discrete momenta appear. For other allowed radii, \B) R is obtained from \B) su , 2 s 
by a similar projection: keep only the components with winding numbers admissible at 
radius R. 

Using ( |2.2tj| ) we can calculate the partition function for other boundary conditions. 
For example, we can replace the Dirichlet state by a Neumann state. In the SU(2) radius 
case we get 

Z S U(2)-{B\q \N)su(2)- 1^ sin(?r | |) ^Xj (q) (2.27) 

i=0,l/2,l,.. V mj V 



and for the infinite radius case we get 

V2 



Z^ 00 = (B\q^ +L ^\N)= ^ 00 (2n\g\)^= X J-( q 2 ) (2.28) 

j=o,i,.. 

where Vq is the SU(2) rotation matrix 

VUM9\) = J { f (1 - J t = sm 2 (n\g\) (2.29) 

We have checked that these formulas agree with the first few orders in their respective 
perturbation expansions. 

A curious result is found if both ends of the string interact with the same sinusoidal 
potential at the self-dual radius. Here the calculation of the partition function is especially 
simple, 

Zj SU{2) — \ B \q \ B ISU{2) — \ iV l e 9 e \^)SU{2) — z 5C/(2) 

Since the SU{2) rotation commutes with the Hamiltonian and annihilates against its ad- 
joint, the partition function is independent of the potential strength and the open string 
energy levels do not feel the potential at all, no matter how strong it is! (This doesn't 
work for other radii where the rotation is operated on by a projection with which it does 
not in general commute.) This is reminiscent of the situation for open strings in constant 
electric fields, where attaching two equal charges to the ends of the string does not shift 
the energy levels [|| . The same feature occurs also in the Kondo model, if one imposes two 
identical Kondo boundary conditions, with a non-zero phase-shift, at the ends of a finite 



strip []18[: The resulting spectrum does not depend on the value of the phase-shift. Such 
phenomena deserve further investigation. 
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3. Calculation and Interpretation of Scattering Amplitudes 

3.1. A Unitarity Paradox and the Role of Solitons 

In this section we analyze the S-matrix for scattering from the dynamical boundary. 
In [12] some exact correlation functions were derived from the constraints of conformal 
invariance. We have just shown that the dynamical boundary state \B) is an SU(2)l 
rotation of the Neumann boundary state \N). We will now show that this remarkable 
result, together with the SU(2)i current algebra of the self-dual c = 1 theory, permits an 
explicit calculation of the entire reflection S-matrix (which turns out to be quite non-trivial 
and to contain some interesting lessons about solitons). 

Let us begin with a description of how one calculates the amplitude for a collection of 
left-moving particles to scatter into a different collection of right-moving particles. We have 
a single field X(z, z) (z = t + ia) defined on the upper half plane. Far from the boundary, 
which runs along real axis, the right(left)-movers are created and destroyed by d z X (d^X), 
and the S-matrix element is the Fourier transform (to pick out the desired in- and out-going 
energies) of (B\d z X(l)d z X(2) . ..d s X(N)\0)/(B\0). If \B) is the trivial Neumann state, 
evaluation is simple: the boundary is eliminated and all the anti-holomorphic fields are 
converted to holomorphic fields by the replacement dX(t, a) — ► dX(t, — a). The resulting 
purely holomorphic matrix element is computed by treating X as a holomorphic free field, 
with propagator (X(z)X(z')) = — log (z — z'), defined on the whole plane. This strategy is 
implicit in the operator Neumann boundary conditions presented in ( |2.15| ) . The distinction 
between in and out fields is now based on whether they are inserted above or below the real 
axis. Now turn on the boundary interactions. The principal result of the previous section 
was that the interaction terms can be rewritten as contour integrals of holomorphic SU(2) 
currents along the boundary (which in fact sum up to a global holomorphic SU{2) rotation). 
The anti-holomorphic fields dX commute with these holomorphic SU(2) currents and 
reflect through the underlying Neumann boundary condition into holomorphic fields in 
the lower half plane just as before. The boundary interactions can be eliminated by 
closing integration contours into the lower half-plane, and carrying out the appropriate 
global SU{2) rotation on every dX operator inserted in that half-plane. For real potential 
strength g, the explicit rotation is 



dX -> cos (2ng)dX + sin(27r#) 



_ e iV2X _|_ e -iV2X 

71 



(3.1) 
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The net result is a very simple prescription: For each incoming field, insert a factor of dX 
above the real axis; for each outgoing field, insert the iS77(2)-rotated version of dX below 
the axis; evaluate the resulting correlator using the free holomorphic propagator for dX 
and Fourier transform to pick out specific incoming and outgoing energies. The nonlin- 
earities introduced by the SU (2) rotation will induce arbitrarily complicated multiparticle 
scattering amplitudes. 

For the 1 — >• 1 correlator this procedure is easily applied. Due to X-momentum 
conservation, only the linear term in (|3.1|) survives to give 

cos 2%g 
(z-z' 



(dX(zrz)BX(z'rz')) = -i^% (3.2) 



Upon Fourier transforming to fix the energy we obtain 

E 

S(E, E') = cos {2itg)-8{E - E') . (3.3) 

The interpretation is that the probability to scatter into only one quantum of the X-field 
is cos 2 (27cg) (the other factors come from our state normalization). We can also use the 
above recipe to evaluate the 1 — > n correlation functions (dX(z)dX(wi) . . . dX(w n )), but 



we find, as already noted in |T2[, that they reduce to a sum of disconnected two-body 
pieces. This means that the probability to scatter into more than one quantum is zero 
and unitarity is violated! We note in passing that a similar paradox arises in the study of 
scattering from a Kondo boundary condition (a fact that has not, to our knowledge, been 
remarked upon). 

The missing probability has presumably leaked into some unaccounted-for sectors of 
the Hilbert space. The classical solution of this system, discussed in Appendix B, shows 
that the classical scattering states indeed lie in disjoint sectors. Initially X(t,0) lies at a 
minimum of the potential, but, after a wave has scattered from the boundary, X(t, 0) may 
be displaced by an integer number of periods to a new minimum. This shift is visible in the 
field X as a topological soliton of winding number n propagating away from the boundary. 
This suggests that the scattering states should be labeled by an integer-valued topological 
charge over and above the usual Fock space variables indexing the number and energy of 
individual particles. The weight-one operators dX(z) we have used to describe asymptotic 
particles are only adequate to describe the trivial topological sector. It seems reasonable 
that we should use the other weight-one operators, e ± l ^ x ( z )^ to generate states in the 
non-zero charge sectors. These exponentials describe solitons (kinks) with a shift of 2it\/2 
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between the values of X far to the right and far to the left of the soliton. Since they are 
holomorphic, they create excitations which travel in one direction at the speed of light and 
can have any energy (just like the excitations created by dX). Unlike the case of massive 
solitons, the profile is in no sense a unique solution of classical equations of motion. 

The well-known arguments about whether degenerate vacuum states survive the tran- 
sition from classical to quantum physics make it less than obvious that such soliton sectors 
are really present in the quantum theory. Soliton sectors are usually associated with de- 
generate vacua. At the classical level, we do have degeneracy associated with the multiple 
minima of the boundary potential. But, the boundary quantum mechanics by itself would, 
of course, have no degenerate vacua. The question is whether coupling it to the external 
field theory allows some memory of the degenerate vacua to remain. But the external field 
theory is that of a massless one-dimensional boson, precisely the kind of theory for which 
infrared fluctuations destroy degenerate vacua. So, the sophisticated quantum mechanic 
would be entitled to conjecture that all memory of degenerate vacua on the boundary 
would be erased by quantum fluctuations and that extra soliton operators would not be 
needed to span the S-matrix. Our explicit construction will show that he is wrong. 

With this in mind, let us look at the scattering of one incident particle into a non-zero 
charge sector. Our rules tell us that the simplest amplitude for scattering into the charge 
±1 sector is 

(dXtfe^*^) - ±-L sin (2^)^-J-^ . (3.4) 

This gives a probability of sin 2 (27T(?)/2 for one charge-zero quantum to scatter into a singly- 
charged soliton state. Obviously, these two single soliton final states by themselves make 
up the missing probability for unitarity. Presumably, the other possible final states, such 
as one soliton plus multiple charge zero states, must have zero probability. In order to 
address such questions, we need a complete orthonormal basis for scattering states that 
spans the soliton sectors. In particular, we need to know whether a soliton-antisoliton state 
is a new object or whether it has already been included in summing over states built out 
of multiple individual charge-zero particles. We will now recast the S-matrix calculation 
in such a way that the complete orthonormal basis for the scattering states, including the 
soliton sectors, is clearly identified. 
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3.2. Infrared Regulators and Orthonormal Bases 

It turns out that there is an infrared subtlety that stands in the way of identifying 
a complete orthonormal basis of states: processes involving solitons are accompanied by 
infinitely many ordinary quanta. To do concrete S-matrix calculations, it seems to be 
necessary to impose an infrared cutoff. Rather than putting the system in a box and 
losing the spatial asymptotic region needed to define the S-matrix, we will impose an 
infrared cutoff by making time periodic with period T. The cut-off restricts all energies to 
be multiples of and only a finite number of states are accessible at any given energy. 
States carry discrete energy labels and, at the end of any calculation, we take the 
continuum limit rii,T — ► oo, keeping Ei = -f-rii fixed. With this cutoff, our base space 
is a half-infinite cylinder. At a — > oo the scattering states are identified with the discrete 
left- and right-moving modes of a "closed string" of length T. The dynamical boundary at 
a = defines an S-matrix by the way in which it reflects left-moving into the right-moving 
states. In what follows we will use our construction of the interacting boundary state to 
evaluate the S-matrix for this discretized problem. 

To simplify the discussion, we consider the case where the field X is compactified at 
the self-dual radius (the ^-matrix at the infinite radius turns out to be identical). The 
holomorphic sector has a level-one SU(2) current algebra, formed by the currents shown 
in (|2.8| ), which can conveniently be used to classify the states. The commutation relations 
are 

3 3^ 3 

+ 3 - - (3 ' 5) 

where = J\ ± J^- Since J 3 (z) = idX(z) / 'y/2, the oscillator part of the Neumann 
boundary state can be rewritten in terms of current operator modes: 

\N) = 2-^exp J2~ J ~nJ-n Yl K-™> • 

n=l m= — oo 

In Section 2, we showed that the dynamical boundary state is generated by a global chiral 
SU(2) rotation acting on \N): 

\B)=e ieaj o\N) , 

where the rotation angle 9 a is determined by the potential strength according to ( |2.23| ). 
This boundary state reflects left-moving modes of the SU(2) current into right-moving 
modes according to 

' J* - J 9 - m ) \B) = (3.6) 
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where (for the case of real potential strength g) 

Jin = e' rj »V!/- irj » = 006(2^) Ji n + ^f^(J_+ n - J_"J (3.7) 



To calculate correlation functions, we use (|3.6| ) and ( |3.7| ) to convert a mixed product of 
left- and right-moving operators into a product of operators of one chirality only. Once we 
have only right- moving operators, we replace \B) by the vacuum, |0), the only component 
of | B) that has no left- moving excitations. Applying this recipe to the 1 — > 1 amplitude in 
the charge zero sector gives 

cos(27rsr) — — = cos(27rsr)-d n _ m . 



(0\B) " »> <0|B> v *'2 

Likewise, one finds that the 1 — > n amplitudes vanish for n > 1. Taking the continuum 
limit, we reproduce our previous calculation, (|3.3|). Since this approach to the calculation 
can be summarized by 

\out) = S\in) , S = e lf)aj Z , 

the S-matrix is a manifestly unitary operator. To verify that in detail for specific examples, 
we have to find an orthonormal basis for the soliton sectors. 

To describe the charge- 1 sector, we may attempt to use the states of the form 

4 J- ni J-n 2 ■ ■ ■ J-nM ■ (3-8) 

Unfortunately, because the J + operator does not commute with the J 3 operators, states 
like this with different numbers of excitations are not orthogonal. The natural basis of 
orthogonal states is instead 

Jl ni J 3 _ n2 ...Jl ni \q) (3.9) 

where \q) is the ground state of the charge-g sector. Since the J 3 commute with each 
other, except for the central charge term, and since the \q) are orthogonal for different q, 
this construct gives a relativistically normalized orthogonal basis for multiparticle states 
in all charge sectors. The \q) are well-defined weight q 2 states in the Hilbert space of the 
SU{2) current algebra. They can be expressed in terms of the action of current algebra 
raising operators on the charge-0 ground state: 

\Q) = J±2 q+ A q+ 3 ■ ■ ■ J±i\0) = e^ X (0)|0) (3.10) 

(note that the sum of the raising operator weights adds up to q 2 as it should). With this 
construction of the basis of states, we can explicitly calculate any S-matrix element. 
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3.3. One-Body Unitarity Check and an Infrared Catastrophe 

Let us now reexamine the scattering of one particle in the uncharged sector from 
the boundary, using the new basis. Following the rules explained above, we find that the 
amplitude for scattering into k particles in the charge- 1 sector is given by 

(0|« Bl ...^ Bfc ^ 1 |0>=i^^* fn . 1 _ En4 , (3-11) 

a matrix element which depends neither on the energies nor on k\ It is easy to see that 
the amplitudes to sectors with \q\ > 1 all vanish. 

The total probability to scatter into a charge one final state can be written as follows 
(taking into account the relativistic normalization of states ( |3.9| )): 

k k=0 ru— nj; 

After writing the delta function as 

»2tt 



the sums over the different 71% decouple. Each sum gives a factor — log(l — e ie ) which is 
then raised to the power k and summed over k to give 



e _21og(l-e ie ) = \_ 



(1-e 



\B\2 



Next we perform the integral over 6 as a contour integral surrounding the origin (avoiding 
the pole at z = 1 by replacing 6 — > 6 + ie in all formulae): 

1 / , z~ m 

® dz- -7T = m . 

2%i J (1-z) 2 



The net result is that the total probability of going into charge one states, ( p.!2|) , is equal to 
sin 2 (27T(7)/2, in agreement with our previous argument and with unitarity. The advantage 
of this more involved calculation is that we are now using an orthogonal basis of states, so 
that there is no doubt whatsoever about proper accounting of probability. 

The absence of energy dependence in the amplitudes ( |3. 1 1|) leads to a kind of "infrared 
catastrophe" for charged final states. This tendency to produce an infinite number of soft 
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quanta is made explicit by rewriting ( |3.11|) in terms of the continuum variables, which 
gives 

<o|^3„, . . . j^JiM ~ i^g&pP - J» . 

Thus, in the limit of T — > oo each amplitude vanishes. Indeed, in the continuum limit it 
would not be consistent to have a nonzero constant amplitude because the probability of 
emitting low energy quanta would diverge (the density of states is ~ -^). However, if we 
introduce a lower cut-off of order 1/Ton the energy and sum over the probabilities, then 
the T-dependence disappears and we are left with a finite result (much like the familiar 
Bloch-Nordsieck calculation in QED) |Tj|]. To show how it works, we sketch the calculation, 
concentrating on the T-dependence and neglecting overall factors. The total probability 
to scatter into k particles is 

E<n*> wi> ~ ^ i % r t r t • • ■ r ^ - ^ ■ ^ 

k k=0 

We decouple the integrals over final energies by using 

S(E - J2 Ei) = ±- / dxe^- E+ T, 

2 71 " J -OO 

and each energy integral now gives 

/ — e lxE = -log(-) + const . 

This logarithm exponentiates and cancels the factor of 1/T 2 that multiplies (|3.13| ) . So, in 
the continuum limit, the probability of emitting any fixed number of particles is zero, but 
the sum over all numbers is finite. Thus, in a sense, the final soliton is built of an infinite 
number of zero energy particles. We should also remark that a physical detector will not be 
sensitive to quanta of energies lower than a certain threshold E m i n , so when we compute 
any scattering cross section we will have to sum over unobserved low energy quanta. As 
that sum will have a similar structure to ( |3.13| ) , the T-dependence will disappear from 



the cross section (the answer will depend, however, on the physical threshold E r 
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3.4- Calculations and Unitarity Checks in Higher Soliton Sectors 



The explicit check of unitarity for a single incoming quantum provides a nice test of 
our new formalism. Another non-trivial test is the 2^2 amplitude in the charge zero 
sector, 

(S)\JlAjlA,\ B ) _ cos 2 (2^) . 
sm 2 (2ng) 

H ^ o mi+m2jni+n2 (mi +TO2 - |ni - - |ni - m 2 \) 

The first term is the disconnected product of two 1 — >• 1 scattering events and the second 
is a connected intrinsic two-body scattering amplitude. Its continuum limit 

— — - — —6 (Ei + Ei — E'i — E' 2 )(Ei + Ei — \E[ — E\\ — \E[ — E2D , 



agrees with the calculation of ref. [12], thus giving a further check on the correctness of 
our formalism. 

We have seen that the expressions for the S-matrix elements depend on the basis chosen 
to represent the states. With a natural orthonormal basis, the amplitudes to produce 
charged states exhibit a classic infrared catastrophe associated with the proliferation of 
soft quanta. Summing the probability over all the states of a given charge, we obtain a 
finite physically meaningful quantity. Below we calculate the total probability to scatter 
to a charge-g final state from an arbitrary multiparticle charge-zero initial state. 

In order to calculate this quantity, it is useful to classify states according to their 
transformation properties under the SU(2) algebra generated by Jq. For instance, any 
one-particle state J*L m \0) transforms as a |l,m) state. The S-matrix is a rotation e l9 J ° 
which commutes with the Casimir operator JqJq- Thus, the S'-matrix conserves the spin 
J and rotates J z according to the spin J representation of e ld J ° . Hence, for a single 
incoming quantum the out-state is also pure spin 1, and the probability to have charge q 
is given by \Vl q (2Trg)\ 2 . This formula agrees with our explicit calculations. 

An iV-particle state transforms as a product of spin-1 representations and we need to 
decompose it into a sum of states with definite spin. Remembering that the original state is 
completely symmetric under particle permutations, we find that it decomposes into a sum 
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2 



of states of spins J = N, N — 2, • • •. For example, a two-particle state |*) = J^_ m J^_ n |0) 
can be shown to decompose as 

I*) = |*2,0> + l*0,0> 

1*2, o) = (^-J- m J- n - -J- m J- n - -J+ m J_ n ^ |o) 

|*0,0> = (^J-mJ-n + \j-mJ-n + \j^ m J-^j |0) 

Using the current algebra, one can show that the probability to be in the |0, 0) state is 

_ 1 4 

3 3m 

(we assumed that m > n) which reduces to 1/3 in the continuum limit (m, n — > oo). The 
probability to be in the |2, 0) state is therefore 2/3 and the probability to find a final state 
of charge q is (in the continuum limit) 

\^, q + 2 -\Vl q (2ng)\\ 

which is totally independent of the energies of the incoming particles. The coefficients of 
the SU(2) representation functions are, as might be expected, the squares of the Clebsch- 
Gordon coefficients for coupling two spin one representations to total spin zero and two. 
For a product of N |1, 0) representations we can write a decomposition 

N 

J 3 -n 1 J 3 -n 2 ---J S -n N \0) = J2 & \j^) 

3=0 

(for even N, only even j are present; for odd N, only odd j are present). The probability 
to find a charge q state as a result of sending in iV charge zero particles is 

N 

J2\C j \ 2 \% q (2ng)\ 2 , (3.14) 

3=0 

where C J is the continuum limit of C J \ and C- 7 is simply the Clebsch-Gordon coefficient. 
In the large n limit, the operators a % n = \J~^\Jn (creation and annihilation operators for 
unit-normalized particle states) have the following commutation relations 



[4, <] = 8^8.. 



m,—n 
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They have vector commutation relations with the SU(2) generators as well. Thus, the 
creation operators a l n (for large n < 0) act as J = 1 c-numbers, and the usual decomposition 
of direct products of c-number tensors gives the Clebsch-Gordon coefficients. 

Eq. ( |3.14f ) clearly shows that an A-particle charge-zero state may scatter only into 
states with charges \q\ < A. An instructive special case is the scattering of a state where 
A particles carry equal energy m, 

\in) = (Jl m ) N \0) , 

into the maximum possible charge q = N. The charge- A component of the corresponding 
\out) state is 

■j^y io> 



sin(27rg) x A 



2i 

and we need to calculate the norm of this state and divide it by the norm of the in-state. 
The algebra is easy to do and we find that 



N-l , , x 

Z 1 \ / 



fc=l 

where C N is the Clebsh-Gordon coefficient for coupling A spin-1 states to total spin A. As 
rii — > oo, C N — > C N as expected. The extra factors in the discretized formula ensure that 
enough energy is available to create a soliton in the charge-A sector: their product vanishes 
for any m such that mN < A 2 , the minimum weight of a charge-A soliton. Although a 
soliton of charge A is a massless particle in the continuum limit, in the discrete case its 
smallest possible weight is A 2 . Thus, it takes a certain minimum energy to create such a 
soliton, but this energy scales to zero in the continuum limit. 



4. Open String Diagonalization by Adsorption Methods 

The discussion in the previous sections dealt with the boundary conformal field theory 
mainly from the closed string point of view. The central object under study was the bound- 
ary state describing the closed string state injected into the worldsheet by the boundary 
interaction. It is also quite instructive to look at this problem from an open string point of 
view. We will show that the boundary interaction can be absorbed into the non-interacting 



open string Hamiltonian by quadrature, giving rise to a " phase-shift" [18[. The method 
is a simple (though somewhat trivial) example of the fusion process that has been used to 
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solve for the multi-channel Kondo fixed point 0. We will use it to derive expressions for 
the partition and correlation functions for comparison with the closed string expressions 
derived in previous sections. In order to avoid certain technical complications, we study 
the case of real coupling g, with the field compactified at the SU (2) radius and a Neumann 
boundary condition at the other end. We leave other interesting cases to Appendix C. 

We start with a free chiral boson X(z), compactified at the SU{2) radius and normal- 
ized so that 

< X(z)X(0) >= -Inz , 
< e PiX{z) e -piX(0) >= Z - P 2 ( 4 ^ 

(we define z = t + ia). The momenta p are restricted by the SU(2) compactification to 
be integer multiples of I/a/2. An antichiral boson X(z) also exists and the two together 
describe a free massless boson on the open line. We confine the system to the line segment 
< a < I by imposing the Neumann condition at both ends. This imposes the following 
relations between chiral and antichiral fields 

X(t, a) = X(t, -a) X(t, l + a)= X(t, I - a) . 

This is equivalent to a single free chiral boson X (z), satisfying periodic boundary conditions 

X(t,a) = X(t,2l + a) 

on a circle of circumference 21. This is just the usual open string. Its partition function 
over a time interval T is 

— 1/24 00 
J v ' Q=-oo 

where w = e~ nT / L and f(w) = YlnLii^ ~ w n ). 

Due to the periodicity of X(z), the set of operators 

j3 (z) = -^d z X(z), J+(z) = e'^W, J-(z) = e~^ x ^ 
V2 

are SU (2)i Kac- Moody generators satisfying periodic boundary conditions. We will use 
them to cast the Hamiltonian for this free open string theory in Sugawara form, the most 
useful starting point for our subsequent discussion of the interacting theory. Because the 
Neumann boundary condition is Virasoro invariant, the chiral and anti-chiral components 
of the stress tensor, T = —^(d z X) 2 and T = —^(dzX) 2 , coincide at the boundaries. This 
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means that the open string hamiltonian density can be written entirely in terms of the left- 
moving fields. By the usual Sugawara arguments for free field theories, it can be rewritten 
equivalently in terms of £7(1) or 5£7(2) Kac-Moody currents: 

H(a) = : J 3 (a)J 3 (a) := : J 1 (a) J 1 (a) := \ : J(a) ■ J(a) : . 

The open string Hamiltonian can thus be written in equivalent ways in terms of the 
Fourier modes J® = (l/2n) f l _ l dae in7r ' 7 / 1 J a (a) of the various currents: 

oo oo oo 



n= — oo n= — oo n=— oo 



It is a straightforward matter to calculate the partition function corresponding to the 
J 3 Sugawara Hamiltonian: For each allowed value, Q, of the charge Jq, there is a 
standard £7(1) Kac-Moody module of conformal weight Q 2 making a contribution of 
u]~ 1 / 24 f(w)~ 1 wQ to the partition function. The allowed values of Q are those of the 
zero-mode momenta of the string and the 577(2) compactification restricts those to be 
integer. Summing over the spectrum of Q gives the standard partition function displayed 
in ( |4.2| ). By symmetry, the spectrum of the other £7(1) charge, Jq, is the same and so is 
the corresponding partition function. 

4-1. Adsorption of the Interaction: Partition Function 



Now we turn on the interaction by adding to (|4.3| ) an interaction, Hi n t, which has 
support only on the boundary a = 0. The main point is that, since the boundary interac- 
tion is critical, the interaction density is a field of boundary dimension one. Since we can 
always recast a conformal open string theory in terms of a single chiral field, the boundary 
interaction must be a chiral field of conformal weight one. The only available such objects 
are the Sugawara currents and the straightforward identification is 

oo 

H int = gj\a = 0)=gj £ J l n . (4.4) 

n= — oo 

For simplicity we take g to be real (in fact it turns out that the phase of g is irrelevant 
for the problem at hand). If we add this to the £7(1) Sugawara Hamiltonian displayed in 
( |1.3| ), we see that the interaction can be absorbed by a c-number shift in the Sugawara 
currents: 

H + H int = j J nJ-n + constant (4.5) 
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Ji = J* + g/2 . (4.6) 

These shift equations, evaluated at n = 0, show that the interaction simply shifts the 
allowed £7(l)-charges Q by g/2, relative to the non-interacting theory. The interacting 
Hamiltonian is the same U(l) Sugawara Hamiltonian as before, except for an overall shift 
of the allowed U(l) charges. The partition function is then a corresponding shift of ( f4.2| ): 

— 1/24 00 

When transformed to closed string variables, this reproduces ( |2.27|) . 

So, as far as the open string energy levels are concerned, the net effect of the interaction 
is a rigid "phase shift" of the spectrum of U(l) charge, which leads to a change in the 
weights of the allowed U(l) Kac-Moody modules. This has been discussed previously in 



the context of the Multi-Channel Kondo model |18[]. It is remarkable that the rich structure 
of the boundary S-matrix which we explored in the previous sections is equivalent to such 
simple open string physics. 

4-2. Adsorption of the Interaction: Boundary Condition 

The simple adsorption method, discussed in the last subsection, can readily be used 
to obtain also the explicit boundary condition, which the SU(2) Kac-Moody generators 
satisfy at the interacting boundary. This is the open string analogue of (|3.6j). 

In order to derive this boundary condition, we note that shifted current operators 
have an analytic time evolution, J x (t, a) = J x (z), due to the shifted form of the Sugawara 
Hamiltonian in ( |4.5| ). This permits us to reintroduce a right- moving current via 

J\t,a) =j\t,-a). 

By definition, the boundary condition on the 1-component of the SU(2) current is 

J 1 {t i 0) = J 1 (t,0). (4.8) 

Since we are only interested in the boundary condition at the a = end of the string, we 
may consider the simpler semi-infinite geometry, by letting / — ► oo. 

In order to find the boundary condition on the other two components of the SU(2) 
currents, we introduce an auxiliary boson field Y satisfying 

J\z) = -^dY(z), J\z) = -^dY(z) , (4.9) 
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and normalized as in ( f4.1| ). Since the relationship between the shifted and unshifted 
current, ( |4.6|) , reads in position space 



J 1 (a) = J 1 (a)+Tvgd(a). (4.10) 

and since l 7 1 (cr) is continuous across the boundary, we see that the auxiliary boson field 
jumps as we cross the interacting boundary: 

Y(0~) - Y(0 + ) = V2ng . (4.11) 

Expressing the remaining components of the Kac-Moody current in terms of auxiliary 
boson field Y, J 3 = =^[e iV2 ~ Y + e~ iV2 ~ Y ], J 2 = ^[e iv ^ y - e~ iV2 ~ Y }, we find immediatedly 
the desired boundary condition: 

f 1 ° ° \ - 

J(t,a)= cos(2tt#) sm(2ng) \ J(t,cr), as a -> (4.12) 

\0 — sm(2irg) cos{2itg) J 
This is the open string analogue of the closed string condition in 



5. Conclusions 

In this paper we have solved the conformal field theory of a free massless scalar 
field living on the half line and interacting via a periodic potential at the origin. An 
577(2) current algebra underlies this system and properties of the interacting theory can 
be computed by 577(2) methods. As the potential strength varies from zero to infinity, the 
theory interpolates between the two trivial limits of a free field subject to Neumann and 
Dirichlet boundary conditions. The intermediate conformal theories are quite nontrivial, 
having complicated multiparticle scattering amplitudes. 

Our methods allow us to compute the S-matrix in complete detail and there are 
some interesting lessons to be learned from the computation. The most important one 
is that scattering states are not labeled just by the number and energies of scattered 
quanta: there is also a discrete charge label, reflecting the freedom of the field at the origin 
to make transitions between the different degenerate minima of the boundary potential. 
These "soliton" sectors of the S-matrix are not visible in naive perturbation theory, but 
must be included in order to maintain unitarity. Other boundary conformal theories, most 
notably those describing the overscreened Kondo model, have a similar unitarity problem 
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which must be resolved, we believe, by scattering into as yet unidentified soliton sectors. A 
further interesting point is that there is an infrared catastrophe associated with scattering 
into soliton sectors: infinite numbers of soft quanta carrying finite amounts of energy are 
always produced. 

These results should have applications in several areas. In the language of open 
string theory, the boundary potential is a spacetime expectation value of the open string 
tachyon field and our results should allow us to construct new open string theory solutions. 
Several condensed matter contexts (Kondo model, quantum wires with impurities) involve 
effectively one-dimensional electrons scattering from localized impurities. Our results will 
describe some critical points of such systems. Finally, the original motivation for this model 
came from dissipative quantum mechanics, and our results provide a complete solution of 
the critical physics of the one-dimensional version of that problem. The simplicity of our 
solution gives some hope that the more interesting critical physics of higher-dimensional 
dissipative quantum mechanics models (for example the dissipative Hofstadter model in 
two dimensions) can be successfully attacked. We hope to return to these matters in future 
publications. 
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Note Added 



Recently, the theory solved in this paper was mapped into a fermionic model [30]. 
All the fermionic calculations agree with the partition functions we have found using the 
boundary state methods. Also, PpJ gives an explicit formula for the R = oo partition func- 



tion of open strings with identical sinusoidal potentials at the two boundaries. This formula 
reveals a "stringy" band structure of energy levels. Below we give a simple calculation of 
this effect using our methods. 

The partition function is given in the closed string channel as Z = (B\q( Lo+L °^\B) R=oo . 
Using the formula for \B) R=ao (2.26), we find 



-ml 



.7=0,1/2,1,.. V m=-j 

where 

V J m ,- rn = (j,m\e 2 ^\ 3 ,-m) 

and we have chosen g to be real. As the strength of the potential varies from zero to 
infinity, g varies from to 1/2. Now, 

E l<-J 2 =ff £ (jMe^ J3 e^e^e-^\jM 

J —IT ^ 7t " 

m=—j m=—j 

sin (2i + 1)0/2 



2tt sin/3/2 
where /3, the angle of the effective rotation, is given by 

sin(/?/4) = cos(7rsr) sin(0/2) . 

We may transform to the open string channel using the formula 

V 1 v fr ,2, sin(2j + l)/3/2 ^ ^ ( „ +&) , 



and we finally find 



2tt /(«> 



n= — oo 



in agreement with [20|. As <fi varies from —it to 7r, varies from — (5 — #) to (| — g). Thus, 
the weights of the open string Virasoro primaries are p 2 — and the allowed values of p 
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form symmetric bands of width 1 — 2g centered on every integer. For g = (no potential) 
the gaps disappear and the spectrum becomes continuous; for g = 1/2 (infinite potential) 
the spectrum becomes discrete, as expected in the tight binding limit. It is remarkable 
that the band structure for open strings is in some ways simpler than the band structure 
for particles. 



Appendix A. The SU(2) Regulator Scheme 

In this Appendix we exhibit the prescriptions for handling ^-function integrals which 
enable us to relate them to the SU(2) algebraic manipulations of section 2. As we showed 
there, in performing the perturbative expansion of Z BD (g,g) we encounter the following 
class of integrals, 



(ei . . . e, 



7 = 1 ^° 



0i(O|2zr) 
6 4 (0\2ir) 



i 



n / n 

k=l J0 l<i<j<n 



Oi(tj 


-U\2iT) 




-ti\2ir) 



(A.l) 



where = ±1. If all e« are equal, the integrals are finite. All other cases suffer from power 
law divergences, and below we construct a natural prescription for handling them. 

From here on we omit the second argument of the theta functions: it will always be 
2ir. The function G(t) = (§jt^j) 2 plays a central role in what follows and we state some of 
its properties: It has a double zero at t = and is periodic under t — > t + 1 and t — > t + 2ir. 
A less obvious, but useful, symmetry property is 



G(t + ir) = 1/G(t) . (A.2) 

Consider now the simplest integral of type (|A.lj ) that needs regularization: (H — ). 
We regularize this integral by shifting the contour of integration away from the real axis, 
t — > t + ie. Once the singularity at t = is avoided, we may shift the contour all the way 
to t — > t + ir (we can shift it freely along the torus because of its periodicity properties, 
as long as we do not cross a pole). 
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FIGURE 1: Shift of the integration contour for the (H — ) case 



Using ( |A.2| ) we find 



(+-) 



gj(0) 
4 (O) 



gUQ) 
4 (o) 



„ U(f) 



(++) (A.3) 



which is finite. We would have obtained the same result if we had shifted the contour by 
— ie, since the residue at t = is zero. 

We will regularize the general integral ( |A.1| ) by shifting the contours of integration by 
an amount determined by their order in the matrix element (£& — > tk + iek): 



(ei ■■■e n ) 



^ 4 (0) 



dti 



dt r 



n 



l<k<l<n 



Oijtk -U+ ie(k - /)) 
9 4 (tk - U + ie(k - I)) 



(A.4) 



FIGURE 2: Shift of integration contours in the regularized n-point function 

In the regulated integral, shown in Fig. 2, the integrands of the individual contour integrals 
are SU(2) current densities, so the A;-th contour should correspond to the action by the 
global SU(2) charge J €k . If this is so, then the regulated integrals ( |A.4| ) are closely related 
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to matrix elements of products of the SU{2) charges, the essence of the result found in 
Section 2. We will show that this expectation is indeed correct. 

First we note that each integral in ( A.4 ) is finite, and has a well-defined limit as e — > 
0+. It value does depend on the order of the £j's, and in calculating the partition function 
it makes sense to average over all possible orderings. It is important to determine how 
the integral is changed by a permutation of two adjacent contours, ti and t i+ i. Consider 
therefore 



(ei... €k H — €k+i-..e n ) 



01(0) 

04(0) 



n+2 



dt\. .dtkdtdt' dtk+i ■ ■ .dt 7 



k n (A.5) 

G(l/-t)- 1 ]lG(t-t i ) e *G(t-t i )- ei J] G(t i -t) ei G(t i -t)-«F(t 1) ..t n ) 

i=k+l 



i=l 



where F contains all the propagators that do not involve the two variables under consid- 
eration (t and t') and tk are complex and of the form tk = Re(tk) +ike. The quantity that 
measures the effect of interchanging two neighboring contours of opposite charge, 



(ei...efc H — efe+i...e n ) — (ei-.-et — Hk+i---^n) 



(A.6) 



has a graphical representation as a contour integral as shown in Fig. 3. 



FIGURE 3: Contour commutation inside an n-point function 



The difference of the two orderings reduces to a contour integral of t' around t, which 
is given by the residue of the integrand in (|A.5| ) at t' = t. Since 



G(t' -t)- 1 



4 (o) y i 
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+ 0((t'-t) ) 



(A.7) 



we need to isolate the term of order t' — t in the remaining factors. After some algebra, 
( |A.6| ) reduces to 



(Ai 



Gfa- t) ** d t G(U-t)-^}F(t 1: .. : t r 

i=k+l 



In performing the integral over t we may use 



f dtG(t + ie)d t G(t + it)' 1 = -2ni 
Jo 



(A.9) 



which finally leads to the result 



(ei...e/s H — efc+x-..e n ) — (ei...ejfc — hefc+i...e n ) = 

'01(0) 



1 5 •••) ^n) 



(A.10) 



71 ^ /t 

= -(2tt) 2 2 Yl e i-2E e « < e i- e ») 

j=fc+l i=l 

This is the formula that will enable us to relate this regularization with the operator results 
of Section 2. 

One may show that the quantities 

(e 1 ...e n )(2 7 rz)- n 

behave under commutation of J £i and J £i+1 in the same fashion as the 577(2) matrix 
elements 



(j,m\J ei ...J €n \j,-m) with m 



En 
i=l C * 



The contour commutation relation ( A.lOj ) corresponds to the SU{2) relation [J + , J~\ = 
2J 3 . This is the essential reason why the perturbative evaluation of Z BD can be reduced to 
the matrix elements of products of the 577(2) charges sandwiched between two boundary 
states. 

Our argument can be made more precise. First we note that there is no doubt about 
the validity of the operator techniques for the case = +1, where the integral is convergent 
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and equal to a matrix element of a power of J + . Now we will relate all integrals to the 
integrals with 6{ = +1. This can be done in two steps. First we reorder the contours to 

bring the integral into the form (+ + ... H ...—). As shown above, the extra terms 

picked up in the process of reordering are equivalent to the commutator terms in the matrix 

element. As a second step, in the integral (+ + ... H . . . — ) we shift the — contours 

up by ir, and using ( |A.2| ) establish the identity (+ + ... H ...—) = (+ + ...+). 




Figure 4: Transformation of (+ + +... ) into (+ + +... + ++) 



The same identity applies to the operator matrix elements, 

aV |(j+)fc(r)r| ^V )=(j '^ |(j+)Hr|j '-^ ) 

Therefore, recursive reduction of the integrals (ei...e n ) to the integrals where all are 
equal, can be fully rephrased in the language of SU{2) matrix elements. This is precisely 
what is behind the remarkable simplicity of the results of section 2. 



Appendix B. Classical Analysis of Scattering 

Our objective here is to gain some intuition about the role of the modes e ±zv ^ x . We 
note that when we scatter dX states, the constant background value of X plays no role. 
Wave packets formed from dX states are localized perturbations that asymptote to the 
same constant value in both directions. That asymptotic constant will of course be at a 
minimum of the boundary potential. However, a localized wave packet can leave X at 
a different minimum of the potential when it reflects from the boundary. The outgoing 
excitation will then be a soliton, interpolating between different asymptotic values of the 
field. Since the degenerate minima of the potential are discrete, the possible soliton charges 
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are integer multiples of some basic unit. The excitations in the soliton sectors are of course 
massless and travel at the speed of light. This is why their quantum description is so 
delicate: locally, they are indistinguishable from the usual zero-charge excitations. 

The equation of motion derived from ( |2.1|) includes the following dynamical boundary 
condition on the field at a = 0: 



1 dX , ^ X(t,0) , . 

(t,0) +27r^sin v ' 1 = (B.l) 



As in the bulk X satisfies the wave equation, we can express it a sum of left an right 
moving waves (with X localized at a particular minimum of the boundary potential): 



X(t, a) = -V2n + f(t + a)+ g(t - a) . 



For an arbitrary incoming wave packet /, we solve for the outgoing wave g by using ( |B.1| ) 
and the condition g(— oo) = 0. We see from (|B.1| ) that X must evolve towards a potential 
minimum for long times, but the final minimum can be different from the initial. 

In figure 5 we display a numerical solution of this problem for a particular choice of 
incoming wave packet such that soliton production actually occurs. This is an example 
of a final state that has to be described by the non-zero charge sectors of the quantum 
SU(2) current algebra. In the classical problem, soliton production either does or does not 
occur for a given incoming wave and it is an interesting problem to identify the incident 
waves for which the outgoing state switches from one value of soliton charge to another. 
We have to be careful not to take the classical calculation too literally: The period of the 
potential is related to the loop expansion parameter of the quantum field theory and the 
particular period for which the theory is at a conformal fixed point corresponds to rather 
strong coupling. 
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FIGURE 5: Wave packet which leaves the field at a different minimum after reflection. 

Appendix C. Open String with Different Boundary Conditions at Both Ends 

In this appendix we treat the case with interactions added at both ends of the open 
string. In particular, we add to the Lagrangian of the non-interacting string, with Neumann 
boundary condition on both ends, an interaction Lagrangian of the form 

/oo 
dt^ ■ J(t, a = 0) + 9 2 -J(t,a = l)}, (C.l) 
-oo 

In order to investigate the effect of the presence of the two boundary conditions, we use a 
conformal mapping. We map the strip < Imz < / ( recall that z = t + ia) conformally 
into the upper half complex ( plane, using 

C = e** (C.2) 

This maps the a = Imz = boundary of the strip onto the positive real £ axis, and the 
a = Imz = I boundary onto the negative real ( axis. 

Due to the transformation law of the the current operators under conformal mappings, 
J(z)dz = j(()d(, this preserves the Neumann boundary conditions. Furthermore, the 
mapping preserves the form of the interaction Lagrangian, now on the real £ axis in the 
half-infinite geometry. However, on the positive real axis we get a boundary condition with 
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interaction parameters 9i, and on the negative real axis one with interaction parameters 

A single boundary condition of this kind has been investigated in Section 4. We 
conclude from ( |4.12| ) that the currents in the half-infinite geometry satisfy 

j(0 = n x j{Q (c.3) 

on the positive real C, axis, and 

M) = 7^7(0 (C4) 

on the negative real ( axis. Here TZi, i = 1,2, represents a rotation with angle 27r|^| about 
the axis 

Both boundary conditions are conformaily invariant and so is the combination of 
boundary conditions described above. Thus we may map conformaily back to the strip. 
This yields 

J{z)=K 1 J{z), lmz = (C.5) 
J(z) =^2 1 J(z), Imz = l (C.6) 

at the two ends of the strip. 

These two equations allow us again to eliminate the right-movers in favor of left- 
movers, defined on a circle of circumference 2/, which satisfy now 

J(a) = HxR^J{<j), as a -> . (C.7) 

This has the same form as the boundary condition in ( [4.12| ), describing a single boundary 
with interaction parameter 6>, which describes the rotation effected by the product IZiR^ 1 
of rotation matrices. The new parameter 9 is easily expressed in terms of 6\ and #2 as 
follows 

-> -> -> -> -> 8 ■ 9 

cos(2tt|6'|/2) = cos (2tt \9\ |/2) cos(27r|6'2|/2) + sin(27r|6' 1 |/2) sin(27r |^ 2 |/2)^- (C.8) 

|0ilr2| 

We have seen in Section 4, that the partition function under the condition ( |(J.7| ) is 
given by (|^) 

Z b iB2 = 1^ w^^ 2 (C.9) 

J \ ' n= — 00 

Notice that we may consider the Dirichlet boundary condition D as a special case of 
_E?2, corresponding to 62 = (1/2,0,0). Thus our general formula, (|C.9| ) , covers also the 
case of the partition function Z BD , for SU{2) radius, chosing 61 = (Re(g), —lm(g),0). 
According to the discussion following (|2.7|) , \9\ = ^ + |, with A given by eq. ( |2.22| ). The 
open string argument just given exactly reproduces the closed string result. 
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